Probing superfluidity of a mesoscopic Tonks-Girardeau gas 
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We study the dynamical response of a Tonks-Girardeau gas on a ring induced by a moving delta- 
barrier potential. An exact solution based on the time-dependent Bose-Fermi mapping allows to 
obtain the particle current, its fluctuations and the drag force acting on the barrier. The exact 
solution is analyzed numerically as well as analytically in the perturbative regime of weak barrier 
strength. In the weak barrier limit the stirring drives the system into a state with net zero current 
for velocities v smaller than v c — nh/mL, with m the atomic mass and L the ring circumference. 
At v approaching v c angular momentum can be transferred to the fluid and a nonzero drag force 
arises. The existence of a velocity threshold for current generation indicates superfluid-like behavior 
of the mesoscopic Tonks-Girardeau gas, different from the non-superfiuid behavior predicted for 
Tonks-Girardeau gas in an infinite tube. 

PACS numbers: 67.85. De,03.75.Kk 
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I. INTRODUCTION 

A typical superfluid displays frictionless flow past an 
obstacle. Superfluid behavior occurs at small flow veloc- 
ities and breaks down at a critical velocity threshold v c , 
where creation of elementary excitations becomes ener- 
getically possible [l|. This phenomenon can be viewed 
in the co-moving reference frame, where impurities or 
tube walls move with respect to the fluid. Following the 
same idea it is possible to probe superfluid-like behav- 
ior by driving a moving barrier with respect to a fluid 
at rest. This "stirring" perturbation has been applied 
to probe superfluidity in a three-dimensional elongated 
Bose-Einstein condensate by moving a blue-detuned laser 
beam inside the condensate 0|. 

For a one-dimensional (ID) fluid, the presence of a 
barrier affects the fluid particularly strongly as particles 
cannot circumvent it; moreover in ID quantum fluctua- 
tions are important. Both aspects are expected to affect 
the transport behavior of a ID fluid induced by a mov- 
ing barrier with velocity v. The possibility of observing 
superfluid-like behavior in ID has been addressed by As- 
trakharchik and Pitaevskii [3j. They have shown that 
for an infinite tube there is no finite critical velocity, but 
there is a non vanishing drag force F& oc v 2K ~ 1 , where K 
is the Luttinger parameter. For a weakly interacting ID 
Bose gas K ^> 1 and the drag force at small velocities 
is very small: the system behaves almost like a 3D su- 
perfluid. In this limit the drag force has been calculated 
using the Bogoliubov approximation 4]. On the other 
hand, in the Tonks-Girardeau (TG) limit of impenetra- 
ble bosons [l,[y| one finds K = 1 (see eg Q) and the drag 
force scales linearly with the velocity as in a normal fluid. 
Calculations of the drag force for a ID Bose gas at large 
interaction strengths are given in [8]. The link between 
drag force and ID superfluidity has been reviewed in [9j . 
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The peculiar transport properties of ID fluids have 
been widely explored in the case of electrons. For an 
infinite fermionic Luttinger liquid the transport behavior 
across a barrier depends dramatically on the strength 
of the Luttinger parameter [lOJ. For K < 1 (eg re- 
pulsive interactions between fermions) the particles are 
perfectly reflected even by a very small barrier - lead- 
ing to an insulating behavior, while for K > 1 they are 
perfectly transmitted even by a very large barrier - dis- 
playing superfluid-like behavior. This change of behav- 
ior as K is increased through K = 1 is also known as 
the Schmid transition [ll(. The case K = 1 (noninter- 
acting fermions) is exactly solvable: for an infinite wire 
the electronic conductance across a large barrier is small 
[lOJ I . hence its behavior is closer to the one of an insu- 
lator. A similar analysis can be done in the bosonic 
case PJi], where the Bose gas with short-range interac- 
tions (K > 1) on an infinite tube is predicted to show a 
supcrfluid-likc behavior at zero temperature, and to un- 
dergo a transition to a non-superfiuid state at K = 1, 
i.e. the TG limit. This transition 11;] can be viewed 
as a localization-delocalization transition of the relative 
phase across the barrier: a superfluid has a well-defined, 
hence localized relative phase, while the insulator has a 
fluctuating, delocalized phase. Superfluidity is destroyed 
by quantum or thermal fluctuations which induce phase 
slips, i.e. jumps of the phase between potential-energy 
minima. 

In the case of a finite ring of circumference L the trans- 
port properties of a ID bosonic fluid - as probed by a 
moving barrier potential - change considerably at small 
velocities. A Luttinger-liquid approach to a weakly in- 
teracting Bose gas [l2[ predicts a superfluid-like behavior 
with critical velocity v c = irh/mL, m being the mass of 
the bosons. The question naturally arises as to what hap- 
pens for strong interactions, especially for K — 1 where, 
according to the analysis of the drag force, the infinite 
system displays nonsuperfluid behavior. 

To answer this question we explore the dynamical 
properties of a mesoscopic TG gas, stirred by a moving 



delta barrier set into motion either abruptly or adiabat- 
ically. This model is of interest since ultracold atomic 
gases trapped in tight ring tr aps are currently being ex- 
perimentally explored [13l - [l9j and the TG gas has been 
experimentally observed in finite ID tubes [20, UM- Us- 
ingan exact mapping onto a noninteracting Fermi gas 
[5ll22J. we can access the exact many-body wavef unction, 
allowing us to analyze the details of its spatial structure 
and the whole dynamical evolution. By calculating the 
integrated current induced in the fluid, the current fluc- 
tuations and the drag force acting on the barrier, as a 
main result we find for the TG gas a superfluid-like be- 
havior with the same critical velocity v c as for the weakly 
interacting Bose gas. 



II. DYNAMICAL EVOLUTION UNDER 
STIRRING DRIVE 

The Hamiltonian of N bosons on a ring of circumfer- 
ence L, stirred by a delta potential U(x, t) = UqS(x — vt) 
is given by 



n 



h b = j: 



2m dx 



—int + ui^^) 



^2g5(xj-xt). (1) 



In the TG regime, corresponding to the impenetrable bo- 
son limit g — > oo, the interaction potential can be re- 
placed by a cusp condition on the many-body wavcfunc- 
tion 



V B (-Xj =x e ...) = 0. 



(2) 



For a gas on a ring we further impose periodic boundary 



conditions ^ B (. 



.) = m B (...Xj + L...) for any j. 



It is possible to obtain an exact analytical solution 
for the many-body wavefunction by mapping the system 
onto a gas of noninteracting fermions subjected to the 
same external (time-dependent) potential [22J, 



% B (xi,...x N ,t) =A(xi,...x N )(l/VNl)det[ipi(x m ,t)} 

(3) 
with A{x\, ...xn) = nj<£sign(izij — xg) being a map- 
ping function which ensures the bosonic symmetry un- 
der exchange of two particles. In each coordinate sector 
x P(i) < xp^2) < ■■■ < £p(7V), where P is a permutation of 
the set {1,2, ...N}, the many-body wavefunction Eq.© 
satisfies the required boundary conditions Eq. (|2J) and is 
the unique solution of the many-body time-dependent 
Schroedinger equation, provided that the orbitals t/)i(x, t) 
satisfy the single-particle Schroedinger equation 



ihd t ipi(x,t) 



2m 



% + Uo6(z-vt))1n(x,t). (4) 




FIG. 1. The energy eigenvalues in units of (irh) 2 /(2mL 2 ) as 
a function of the stirring wavevector q in units of n/L. Each 
parabola represents a state of different angular momentum. 
These states are coupled by the barrier that opens gaps at 
q n — nn/L and thus allows to change the branch of angular 
momentum. The inset shows the first avoided level crossing 
between the states with zero and one quantum of angular 
momentum. 



choose as the initial condition the ground state orbitals of 
a TG gas in the presence of a nonmoving barrier localized 
at x = H, i.e. in Eq.© we set ipi(x,0) = ^(x), 
where for spatially even orbitals we have 

^ 0) (x) = ^cos{k^(\x\-L/2)] (5) 



a/-; 



with wavevectors fixed by the transcendental equation 
k\ tan(fc ; L/2) = mUo/ti 2 and normalization given by 

A/"/ 0) = ^/2L(l + sin(fc / (0) L)/fc, (0) L), while for odd or- 
bitals we have 



(j)\°\x) =i(-l)\j^rsia[2Trlx/L], (6) 

with I — 1,2,3... integer, and the phase factor in Eq.© 
has been chosen for consistency with the forthcoming 
Eq.CJJ. 

The solution for the dynamical evolution according to 
Eq.© then reads 

Mx, t) = e^e"^ 2 '/ 2 ™ J2 Cjie- lE ^ /h ^, q {x - vt), (7) 



where hq — mv is the barrier momentum. The time- 
independent overlaps Cji contain the information about 
the initial condition and arc defined as 



Cjl = <^ 9 |e-^|^ (0) > 



(8) 



Exact Solution 



The orbitals <j)j tq (x) and the energies Ej jQ are obtained 
from the stationary Schroedinger equation in the frame 
co-moving with the barrier, 



If the barrier is set into motion instantaneously at t = 
we have found an analytical solution for Eq.© [23|]. We 



-—di + U S(x) ) cf> M = E hq ^ q (9) 



with twisted boundary conditions (TBCs) 4>j. q (x + L) = 
e~ lqL <j)A q {x). The solution of Eq.© is obtained by ex- 
pansion with respect to plane waves, 



<t>j,q(x) 



The 
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iq-2 ( e i{kj{x--2) ^ £ e -ikj(x--z)\£ Jq^ 



normalization 



factor 



Mi 



(10) 



L{\ + A 2 + 2Aj sin(kjL)fkjL) and the amplitude 

Aj = sin[(fcj + q)L/2]/ sin[(kj — q)L/2] are functions 
of the wavevectors kj given by the solution of the 
transcendental equation 



mUo sm(kjL) 

h 2 cos(qL) — cos(kjL) 



(11) 



This also fixes the energy eigenvalues Ej = h 2 k 2 /2m that 
are shown in Fig|TJ Finally, using Eq. ([TT))) and respec- 
tively Eqs.© and © above, it is possible to find an 
analytical expression for the overlap coefficients, which 
read 
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for the overlap with an initial even orbital, and 
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for the overlap with an initial odd orbital, with Jq(x) 
sin(x)/a; being the first spherical Bessel function. 
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FIG. 2. (Color online) a) Time-averaged, integrated particle 
current in units of (TvNh)/(mL 2 ) and b) current fluctuations 
in units of (Nhkir)/(mL), kF being the initial wavevector 
of the Nth particle, for a TG gas (red dotted line) and a 
noninteracting Bose gas (cyan dashed line) subjected to the 
nonadiabatic stirring described in the text, as a function of the 
stirring wavevector q in units of n/L, for a barrier strength 
rnUoL/h 2 — 1 and N — 3 particles. The behavior of the 
TG current following an adiabatic switching on is also shown 
(black solid line). Insets: Zoom on the second current (current 
fluctuation) peak. The numerical results (same notation and 
color codes as the main panel) are compared with the results 
of the perturbative approach (thin solid lines). 



B. Perturbative solution for weak barriers 

In order to develop a qualitative understanding of the 
behavior of the system under stirring, we study the case 
of weak barriers, i.e. for XL < 1, where it is possible 
to obtain the dynamical evolution of the wavefunction 
([7]) perturbatively. Due to the ring periodicity the prob- 
lem of the time- independent Hamiltonian in Eq. © with 
solutions that obey TBCs is equivalent to the one for 
a particle in a periodic potential and we can therefore 
expand its solutions in Bloch waves 



4>j, q (x) 



VI 



P6Z 



27rpx/ L 



.,W) 



(14) 



with eigenvalues Ej, q and coefficients o ' ! 2 ir /l> •?' ^ e " 
ing the band index. Indeed the eigenstates <j>j^ q (x) of 
the Hamiltonian in Eq.© are written as a product of a 
plane wave exp(— iqx) and a lattice- (in our case ring-) 
periodic function and thus obey Bloch's theorem. The 



parameter q plays the role of the quasi-momentum. It 
provides a periodic continuous parametric dependence of 
the states and the energy on the stirring velocity. If the 
barrier is weak it will only couple two states 2ira/ L and 
27T&/L, a, b G Z of different angular momentum depend- 
ing on the concerned avoided level-crossing, see FigUJ 
Consequently, at fixed avoided level crossings around 
q n = rnr/L, n being an odd integer only two elements 
of the sum in Eq. (fhll) contribute and it reduces to 
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for odd values of j. The expansion coefficients a _ 2 , L 
are readily obtained by solution of the effective two-level 



problem, 
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(17) 
(18) 

(19) 



the +(-) solution referring to u q (v q ), and we have in- 
troduced 5q = q — q n , the deviation of q from q n . Equa- 
tions (fl5l) and (fl6]) also hold for even values of n choosing 
j even |25J . 

We finally obtain the full wavefunction at the avoided 
level crossing inserting Eq. (fl"5"j) and (fT6| in Eq.(|7|), 



e iqx e -iq 2 t/2m 
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(20) 



with 



i E 3\q> E j+i,g} 
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(21) 



being the energy of the jth and the (j + l)th band. Here 
the Yl i s over the even j or the odd j only, depending on 
n. The coefficients Cji are determined using Eq. (|S|) where 

the initial orbitals are plane waves </>J = exp(ikix) / \/L 
with ki = —tv (I — 1)/L for odd I and ki = irl/L for even 
I, yielding 



Cji = v q J 'dj. n+ [-i 
- „<J) 



u q""j,n+l-l 



Cji 



u a u t.n— 



CO; 



j,n-l 



c j+lJ — u q &j,n-l 



for I odd 



for I even 



(22) 



III. INTEGRATED PARTICLE CURRENT 

The knowledge of the many-body wavefunction allows 
to derive several observables of interest. We focus first 
on the integrated particle current / = (1/ ' L) J dxj{x,t), 
which for the ring geometry is proportional to the angular 
momentum. The current-density j (x, t) coincides with 
the one of the mapped Fermi gas [26| , 



iV 



j(x,t) = — Im.y^jp*{x,t)d x il>i(x,t). 



(23) 



Using this expression we evaluate the integrated current 



1=1 



dq \ 2 



+ Im£«4q,e^-^% ,(24) 

i& J 

where the amplitudes Fij — (<f>i\d x \4>j) are given by 



F„ = ~ tk : k i L2 JoKk - k 3 )Ll2]J Q {{ki + k 3 )L/2] 



jViMj 



sixi{qL) 



sin[(fc l - q)L/2] wa[{kj - q)L/2] 



(25) 



The result (IMt for the spatially integrated current dis- 
plays a time-independent contribution which has the 
form of a generalized thermodynamic relation / ex dE/dq 
and a time-dependent oscillating term, with several fre- 
quency components originating from the multimode na- 
ture of the state of the TG gas. In order to gain further 
insight we focus on the case of weak barrier strength. 
Close to any level crossing, Eq. (|20| allows to obtain an 
explicit expression for the integrated particle current, 



A' 



I(t) = —J2^n~k l ) Uq ^V q ^ 



1=1 



where we introduced 



1 - cos(2AE^'H)} , (26) 



Jnl — 



\n-l\ 
(n + l-1) 



for 
for 



I even 
J odd 



(27) 



Here AE, 



(0 



E, 



l+l,q 



Ei iQ is the energy splitting of 



/=i 



the effective two-level system that governs the time- 
dependence of the current. If averaged over sufficiently 
long times T > K/ AE q , this dependence vanishes. In 
FigGk we illustrate the resulting time-averaged current 
for a weak barrier. We find that it displays a narrow max- 
imum for values of the stirring momentum hq n that are 
equal to integer multiples of hq c — hir/L corresponding to 
a critical stirring velocity v c — nh/mL. This result shows 
that it is difficult to stir the mesoscopic TG gas except 
at velocities that are multiples of v c . At generic values of 
the stirring momentum away from hq n the fluid belongs 
to the zero current branch, which for nonadiabatic stir- 
ring momenta hq > Hq c corresponds to an excited state. 
The occupation of such highly excited states is due to 
the instantaneous switching on of the barrier motion. At 
the special values of the stirring momentum hq n avoided 
crossings between single-particle branches with different 
angular momentum occur. This allows for the occur- 
rence of superpositions of states with different angular 
momentum [23j which yield a non zero average current. 
In contrast to the nonadiabatic stirring mechanism, an 



adiabatic switching on leaves the system in the ground 
manifold {fey} with j = 1..JV and corresponding energy 
E g (q) = Ylj=i fo 2 k 2 /2m. In this case the integrated, time 
averaged current is given by the thermodynamic relation 
/ = (Nhq/m — dEg/dhq)/ L and has a staircase behavior 
as a function of the stirring momentum, also illustrated 
in FigHh, 

Figure [2^, shows for reference the corresponding time 
averaged current of an ideal Bose gas subjected to the 
same stirring drive. For momenta close to Hq c we find 
that the current induced in an ideal Bose gas is larger 
than the one of the strongly interacting TG gas. This ap- 
parently surprising result is readily explained by noticing 
that the ideal Bose gas occupies the lowest single-particle 
level which is the most affected by a weak delta-barrier 
potential since u q and v q show a very smooth depen- 
dence on q. For large stirring momenta hq 3> hq c this ef- 
fect vanishes since then also the ideal Bose gas is in a very 
highly excited state. In the inset of Figj2^, we show that 
a very good agreement is found at weak barrier strength 
between perturbation theory and exact calculation for 
momenta close to hq n . 

Finally, we remark that our results for the integrated 
current are also in agreement with the prediction of non- 
superfluid behavior of a TG gas in an infinitely long wire. 
In that case, L tends to infinity keeping N/L constant, 
and the critical velocity v c tends to zero. The observed 
supcrfluid effect is truly mesoscopic, i.e. associated with 
the finite size of the ring. 



IV. PARTICLE CURRENT FLUCTUATIONS 

To further explore superfluid behavior, we complement 
the study of the integrated particle current by analyzing 
its fluctuations. This quantity measures how well de- 
fined the particle current is with respect to unavoidable 
quantum fluctuations. The time-averaged particle cur- 
rent fluctuations are defined as 



AI = V(I 2 )-(I) 2 , 



(28) 



where the time-averaged integrated particle current is 
given by (I) = l/T J dtl(t). The second moment of 
particle current reads 

(/ 2 ) = tV^ II dxdy C 2 (x,y,t) , (29) 

Ll " 1 J-T/2 J-L/2 




b) 



(NO 










- 




/P^ 


\ 


\£y 




\*v 






- 



2 4 6 



t/to 




FIG. 3. (Color online) The drag force exerted by the fluid 
on the moving barrier in units of Uopo, po — N/L being the 
average density of the fluid, as a function of time in units 
of to = mL 2 /ivh for a TG gas (red dotted curves) and an 
ideal Bose gas (cyan dashed curves) for three values of stirring 
wavevector, q = 0.25 q c , (a), q = 0.925 q c (b) and q = q c (c), 
for N — 3 particles and a barrier strength mlloL/h 2 — 1. 
The numerical results are compared with the results of the 
perturbative approach (thin solid lines). In panel (c) for q = 
q c the analytical curves for ideal bosons and the TG gas are 

identical since 5q = and therefore AEq and <u q ,Vq > are 
independent of I. 



where 



C 2 {x,y,t) = {^B\j{x,t)j{y,t)\^ 1 



(30) 



is the spatially integrated equal-times current-current 
correlation function. 

The particle current fluctuations derived from pertur- 



bation theory read 

2 N 
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(31) 



FigEb shows the time-averaged particle current fluctua- 
tions for a non-adiabatic stirring. For stirring momenta 
far away from hq n the qualitative behavior of the fluc- 
tuations is similar to the one of the current itself. Con- 
sequently in the superfluid regime, away from q n there 
is neither an induced current nor strong current fluctua- 
tions. The quantum state of the system corresponds to a 
well-defined angular momentum state which explains the 
smallness of the fluctuations. On the other hand, around 
the current peaks, i.e. for stirring momenta hq n , we ob- 
serve strong fluctuations. At these special points, the 
system oscillates between two equally possible angular 
momentum states, leading to the appearance of strong 
current fluctuations. At q n we further observe the for- 
mation of a dip in the fluctuations, which can be under- 
stood qualitatively as follows. The nonadiabatic current 
in perturbation theory is given by the weight combination 



4°- 



,(0 



Ua'Vq' , a structure similar to the second cumulant char- 



acteristic for the partition noise for particles tunneling 
across a barrier [27] . Its fluctuations are consequently 
given by the fourth cumulant of the same distribution 
which gives these splittcd peak structures. 



V. DRAG FORCE 

The drag force |3J acting on a barrier potential Vbarr{x) 
is defined as Fd — ~(d x Vb ar r(x)) where (...) denotes the 
quantum average with respect to the state of the system. 
For the specific case of a delta barrier potential Uq8(x) 
the drag force is simply related to the particle density 
profile n(x), according to 

Fd = U -[d x n(x)\ x=0 + + d x n(x)\ x=0 -}. (32) 

For a TG gas the density profile is readily calculated 
within our exactly solvable model, as it coincides with 
the one of a noninteracting Fermi gas, 



JV 



n{x,t)=J2\Mx,t)\'' 



(33) 



The expression for the drag force acting on a delta barrier 
potential suddenly set into motion across a TG gas in 
perturbation theory reads 



F, 



L ^ 

i=i 



(qa-k) (utf»«)t;tf"»>) sin(2A£< /B,) t) • 



(34) 



Fig[3] shows the drag force obtained from both the ex- 
act model Eqs. ([32]) . ([33]) and and from perturbation 
expansion at weak barrier strength. It is an oscillating 
function of time with a typical frequency associated with 
the energy-eigenvalues separations AEq whose main 
contribution comes from the Fermi level, E^+i, q — EN,q- 
Notice that since the energy levels depend on the stirring 
momentum the typical oscillation frequency is also mo- 
mentum dependent. At increasing stirring momentum, 
we find that the behavior of the drag force closely follows 
the results for the integrated current presented in the pre- 
vious section. In the case of hq < hq c the drag force is 
vanishingly small, confirming the picture of a superfluid 
mesoscopic TG gas. If the barrier momentum is instead 
chosen close or equal to Hq c , then it is possible to trans- 
fer angular momentum to the fluid and the drag force 
is nonvanishing, illustrating the breakdown of superflu- 
idity. We find the same velocity threshold as the one 
predicted in Ref. 12] using a Luttingcr-liquid description 
of a weakly interacting Bose gas. Finally, for momenta 
larger than Hq c except at the special momenta close to 
integer multiples of Hq c , we find that the drag force van- 
ishes as the system is in a very peculiar, highly excited 
state which does not couple to the barrier motion. 

In Fig|3]we also show the drag force originating from 
the stirring of an ideal Bose gas as a function of time. 
The typical oscillation frequency of the drag force of an 
ideal Bose gas is fixed by the first energy level splitting 
E2, q — E\ t q. Due to the flatter dispersion of the first en- 
ergy levels as compared to those at the Fermi surface, for 
the ideal Bose gas the dependence of the drag-force os- 
cillation frequency on the velocity is weaker than the one 
found for the TG gas. The three panels of Fig|3] corre- 
spond to increasing values of the stirring momentum. At 
small momenta also the ideal Bose gas shows mesoscopic 
superfluidity. For momenta close to Hq c the drag force of 
the ideal Bose gas is larger in magnitude than the one of 
the TG gas, illustrating the more important effect of the 
barrier on the ideal Bose gas. For a stirring momentum 
equal to hq c the drag force of the TG gas and of the ideal 
Bose gas are equal in magnitude, as it is the case of the 
integrated current (see again Fig[2^). 



VI. CONCLUSIONS 

In summary, using an exactly solvable model, we have 
investigated the transport properties of a TG gas on a 
ring trap stirred by a localized barrier potential moving 
along the ring circumference. The calculation of the spa- 
tially integrated particle current its fluctuations and of 
the drag force acting on the barrier shows the existence 
of a velocity threshold for setting the TG gas into mo- 
tion, at a value given by v c — hir/mL. The same value 
has been predicted for the weakly interacting mesoscopic 
Bose gas by a Luttinger- liquid calculation [12J. We have 
also shown that depending on the way the final velocity 
of the barrier is reached, different quantum states are ac- 



cessed. In particular, we have found that if the barrier 
motion is suddenly switched on, it is difficult to trans- 
fer angular momentum to the system except close to v c 
or its integer multiples. Finally, by comparing with the 
corresponding results for a noninteracting ID Bose gas, 
we have shown that the stirred current, the current fluc- 
tuations and the drag force depend on the interaction 
strength. These limiting cases allow us to obtain a qual- 
itative picture for the behavior of Bose gases for both 
very weak and very strong interactions. Given the ex- 
perimental progress in recent years, the realization of a 
experimental set up that allows to study the mesoscopic 



TG gas on a ring seems close to reach. In perspective, it 
would be interesting to develop tools for describing the 
case of a ID Bose gas with arbitrary interactions sub- 
jected to an out of equilibrium drive. 



ACKNOWLEDGMENTS 

We acknowledge discussions with Denis Basko, 
Edouard Boulat, Roberta Citro and Helene Perrin. We 
thank the PEPS-PTI project "Quantum gases and con- 
densed matter", the MIDAS STREP project and the 
Handy-Q ERC project for financial support. 



[1] E. M. Lifshitz and L. P. Pitaevskii, Statistical Physics 

Part 2 (Butterworth-Heinemann, Oxford, 1980). 
[2] C. Raman et al., Phys. Rev. Lett. 83, 2502 (1999). 
[3] G. E. Astrakharchik and L. P. Pitaevskii, Phys. Rev. A 

70, 013608 (2004). 
[4] A. G. Sykes, M. J. Davis, and D. C. Roberts, Phys. Rev. 

Lett. 103, 085302 (2009). 
[5] M. D. Girardeau, J. Math. Phys. 1, 516 (1960). 
[6] M. D. Girardeau, Phys. Rev. 139, B500 (1965). 
[7] M. A. Cazalilla, J. Phys. B 37, SI (2003). 
[8] A. Y. Cherny, J.-S. Caux, and J. Brand, Phys. Rev. A 

80, 043604 (2009). 
[9] A. Y. Cherny, J.-S. Caux, and J. Brand, preprint 
larXiv:1106.6329l 
[10] C. L. Kane and M. P. A. Fisher, Phys. Rev. Lett. 68, 

1220 (1992). 
[11] A. Schmid, Phys. Rev. Lett. 51, 1506 (1983). 
[12] H. P. Buchler, V.B. Geshkenbein and G Blatter, Phys. 

Rev. Lett. 87, 100403 (2001). 
[13] S. Gupta et al., Phys. Rev. Lett. 95, 143201 (2005). 
[14] A. S. Arnold, G S. Garvie, and E. Riis, Phys. Rev. A 

73, 041606 (2006). 
[15] O. Morizot et al, Phys. Rev. A 74, 023617 (2006). 



[16] K. Henderson, G Ryu, C MacCormick, and M. G 

Boshier, New J. Phys. 11, 043030 (2009). 
[17] A. Ramanathan et al., Phys. Rev. Lett. 106, 130401 

(2011). 
[18] B. E. Sherlock et al, Phys. Rev. A 83, 043408 (2011). 
[19] S. Moulder et aL. larXiv:11 12.0334v2 (2012). 
[20] T. Kinoshita, T. Wenger, and D. S. Weiss, Science 305, 

5687 (2004). 
[21] B. Paredes et al, Nature 429, 277 (2004). 
[22] M. D. Girardeau and E. M. Wright, Phys. Rev. Lett. 84, 

5691 (2000). 
[23] C Schenke, A. Minguzzi, and F. W. J. Hekking, Phys. 

Rev. A 84, 053636 (2011). 
[24] M. D. Girardeau and A. Minguzzi, Phys. Rev. A 79, 

033610 (2009). 
[25] For even n the first band j = 1 is in a definite state of 

angular momentum and its wavefunction reads <f> q (x) = 

exp(— iqx) exp(iTrnx / L) / \[L. 
[26] R. Citro, A. Minguzzi, and F. W. J. Hekking, Phys. Rev. 

B 79, 172505 (2009). 
[27] Y. M. Blanter and M. Biittiker, Phys. Rep. 336 1, (2000). 



